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Abstract: Results are presented from radioactivity screening of two models of photo- 
multiplier tubes designed for use in current and future liquid xenon experiments. The 
Hamamatsu 5.6 cm diameter R8778 PMT, used in the LUX dark matter experiment, has 
yielded a positive detection of four common radioactive isotopes: 238 U, 232 Th, 40 K, and 
60 Co. Screening of LUX materials has rendered backgrounds from other detector mate- 
rials subdominant to the R8778 contribution. A prototype Hamamatsu 7.6 cm diameter 
R11410 MOD PMT has also been screened, with benchmark isotope counts measured at 
<0.4 238 U / <0.3 232 Th / <8.3 40 K / 2.0 ±0.2 6() Co mBq/PMT. This represents a large 
reduction, equal to a change of x ^ 238 U / x \ 232 Th / x \ 40 K per PMT, between R8778 
and R11410 MOD, concurrent with a doubling of the photocathode surface area (4.5 cm 
to 6.4 cm diameter). Co measurements are comparable between the PMTs, but can 
be significantly reduced in future R11410 MOD units through further material selection. 
Assuming PMT activity equal to the measured 90% upper limits, Monte Carlo estimates 
indicate that replacement of R8778 PMTs with R11410 MOD PMTs will change LUX 
PMT electron recoil background contributions by a factor of after further material 
selection for 60 Co reduction, and nuclear recoil backgrounds by a factor of x ^. The strong 
reduction in backgrounds below the measured R8778 levels makes the R11410 MOD a very 
competitive technology for use in large-scale liquid xenon detectors. 

Keywords: Detector modeling and simulations I, Photon detectors for UV, visible and 
IR photons (vacuum). 
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1. Introduction 

Photomultiplier tubes (PMTs) are commonly used for detection of photon signals in scin- 
tillation detectors. In the case of low-background detectors, the proximity of the PMTs to 
the active scintillation region can lead to a background dominated by PMT radioactivity. 
Understanding PMT radioactivity and selecting PMTs based on radiopure material com- 
position is therefore of primary importance for experiments making use of low-background 
detectors, as the use of ultra-low background PMTs can greatly extend background- free 
running times. 

The LUX dark matter direct detection experiment uses 350 kg of liquid xenon for the 
detection of recoils resulting from the scattering of Weakly Interacting Massive Particles 
(WIMPs) U[. The experiment uses 122 Hamamatsu R8778 PMTs directly above and below 
the active region to detect scintillation light from xenon nuclear recoils [0. R8778 PMTs 
are well suited for use in liquid xenon detectors. Xenon scintillation light is narrowly peaked 
at 178 nm [0, @J; the LUX R8778 PMTs feature a measured average 33% quantum efficiency 
and 90% collection efficiency at this wavelength. They offer single-photon sensitivity with 
a resolution of 30%. Their 4.5 cm diameter photocathode allows for high surface area 
coverage per tube compared to PMTs used in previous xenon dark matter detectors, such 
as the 2.2 cm R8520 |5|, 0, 0] . The R8778 was also specifically designed to operate at liquid 
xenon temperatures, in the range 165-180 K. Extensive testing has been performed on the 
LUX PMTs, and their behavior in liquid xenon is well understood. 
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Through extensive screening and material selection, it has been possible to render the 
background contributions from other LUX construction materials subdominant to that of 
the PMTs UK]. The PMTs contribute backgrounds in the form of gamma rays and neutrons 
produced by radioactive decay, as well as neutrons from (ct,n) reactions in PMT materials. 
Neutron backgrounds are of particular importance to dark matter detectors due to their 
similarity to a WIMP nuclear recoil signal. Precise measurements of all radioactive isotopes 
must be made to ascertain the impact of the PMTs on LUX backgrounds. 

LUX has initiated a development program in conjunction with Hamamatsu to identify 
candidate PMTs for use in future dark matter detectors. The LUX-ZEPLIN (LZ) liquid 
xenon detector will be a factor of x20 larger in mass than LUX, with corresponding surface 
area increase of x24 ||. The PMTs used in these detectors will ideally be larger than the 
R8778, and offer lower radioactive backgrounds per unit area. 

The development program has resulted in the creation of several large-area photo- 
cathode PMTs, notably the 7.6 cm diameter R11065 and its low-background variant, the 
R11410 MOD @, HO). The R11065 was tested at Brown University and was found to yield 
an identical response to xenon scintillation light as the R8778, with gains measured up 
to 10 7 at 1500 V and a single photoelectron resolution of 32% while immersed in liquid 
xenon at 175 K. The R11410 MOD PMT has been selected as a candidate for use in the 
LZ program. The R11410 MOD and R8778 PMTs are shown for comparison in figure [I]. 

2. Counting Results 
2.1 SOLO 

The Soudan Low-Background Counting Facility (SOLO) is located at the Soudan Under- 
ground Laboratory, at a depth of 2.0 km.w.e JTTj]. The SOLO geometry is shown in figure 0. 
The counting chamber uses the "Diode-M" detector, a 0.6 kg high-purity germanium de- 
tector housed in a 0.15 cm thick copper shield. The 20 cm x 25 cm x 30 cm chamber is 
enclosed with a minimum of 30 cm lead shielding on every side. The inner 5 cm lining of 
the chamber is comprised of ancient lead, with 210 Pb activity measured below 50 mBq/kg. 
A mylar shell and 2.5 slpm nitrogen gas purge are used to eliminate gaseous radon from 
the chamber. Detector output is shaped with a 2 jUs shaping time, using a Spectroscopy 
Amplifier, before digitization at a dedicated acquisition computer. Energy bin counts are 
written to data files every four hours for analysis. 

SOLO has been used since 2003 for screening and cataloging of low-background ma- 
terials for use in liquid xenon detectors. SOLO was used in the screening program for 
XENON10 construction materials, and has been heavily engaged in the LUX material 
screening program since 2007 JH, 0, O, E3 • 

Detector sensitivity calibrations are accomplished using Geant4 Monte Carlo simula- 
tions. The geometry of the entire SOLO detector and chamber are reproduced, as well 
as the geometry and location of the counted sample. Emission spectra for each screened 
isotope, including the full decay chains of 238 U and 232 Th, are reproduced. These simula- 
tions account for effects including energy-dependent detection efficiency, sample location, 
sample size, and sample self-shielding. The simulations produce scaling factors between 
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Figure 1. The Hamamatsu R8778 PMT, used in the LUX experiment (left), and the R11410 MOD 
PMT, a candidate for use in the LZ detector (right). The R11410 MOD features a 6.4 cm diameter 
photocathode, offering twice the surface area of the 4.5 cm diameter R8778 photocathodc. 

detector counts and source emission intensity for each gamma line, referred to as the de- 
tection efficiency £. Analytic estimates of £ can be generated to cross-check simulation 
results. However, it is very difficult to include the effects of finite sample size and sample 
self-shielding in these estimates, whereas those effects are automatically included in the 
Monte Carlo analysis. 

A typical background spectrum with the chamber empty is shown in figure |3| (black) 
for 21 live days. Features are seen at 662 keV and 1461 keV, corresponding to gamma 
lines from 137 Cs and 40 K, respectively, and are used to verify energy bin calibration in 
conjunction with the 46.5 keV line from Pb decay. Background baseline rise below 
1 MeV is caused primarily by Bremsstrahlung radiation from 210 Pb beta decay in the inner 
lead shielding layers. Initial detector calibrations were performed using several samples 
with known contamination. Backgrounds are remeasured periodically to check for chamber 
contamination or radon leakage, and are used to periodically verify calibrations. 

2.2 PMT Activity Measurements 
2.2.1 R8778 

A total of 20 R8778 PMTs were screened at SOLO in order to yield high-statistics mea- 
surements of their radioactivity. Ten R8778 PMTs, purchased for initial testing at Brown 
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Figure 2. The open SOLO chamber. The 0.6 kg high-purity germanium detector is housed in a 
1.5 mm thick copper shield. The chamber is lined with >30 cm of lead shielding. The inner 5 cm 
lead layer was selected for its measured 210 Pb content (<50 mBq/kg). 

University, were screened in two separate batches of five. After the purchase of the 122 
PMTs for LUX, two additional batches of five LUX PMTs were counted. Average results 
are listed in table 0, with the screening spectrum from one of the batches of five R8778 
PMTs shown in figure |j| Results indicate the presence of benchmark isotopes 238 U, 232 Th, 
40 K and 60 Co in amounts of 9.5 ±0.6, 2.7 ±0.3, 66 ±6, and 2.6 ±0.2 mBq/PMT, respec- 
tively. The measured values represent an overall gamma ray emission rate per PMT that 
is a factor of 3 below the conservative estimate originally used in the LUX proposal, which 
used the manufacturer-supplied upper limits of <18 238 U / <17 232 Th / <30 40 K / <8 60 Co 
mBq/PMT. 

2.2.2 R11410 MOD 

An R11410 MOD PMT mechanical sample was screened at SOLO for a period of 19 
live days. This PMT was produced after a thorough material screening program by the 
manufacturer, with strong reductions in benchmark isotopes expected. Placement of the 
R11410 MOD was optimized within the chamber, and the counting time was extended a 
factor x2 beyond the typical time for the R8778 batches, in order to maximize sensitivity 
for the single sample. 

The R11410 MOD showed a very significant reduction in 238 U, 232 Th and 40 K content 
compared to R8778 counting results, and the signal for these isotopes is consistent with 
background. The only positive signal was found at 1173 and 1333 keV, consistent with the 
presence of 60 Co in amounts comparable to those measured in the R8778 PMTs. Further 
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Table 1. R8778 and R11410 MOD radioactivity screening results, obtained at SOLO. All PMTs 
were screened for four benchmark isotopes identified by the factory: 238 U, 232 Th, 40 K, and 60 Co. 
R8778 PMTs were counted in four batches of 5 units, with an average 9 live days per batch. 
Averaged results are listed. The R11410 MOD sample was counted for 19 live days. Separate 
columns are provided for completeness for "early" and "late" 238 U chain measurements, allowing for 
potential equilibrium breakage primarily due to 226 Ra solubility in water; however, it should be 
noted that there is no a priori reason to assume an equilibrium breakage. Errors are statistical and 
quoted at ±lo~. Upper limits are given at 90% confidence level. Results are normalized per PMT 
for comparison. 



PMT 




Activity 


mBq/PMT] 








238-lj ( 234m Pa) 


238 U ( 226 Ra) 


232 Th (228 Ra) 


40 K 


60 Co 


R8778 


<22 


9.5 ±0.6 


2.7±0.3 


66±6 


2.6 ±0.2 


R11410 MOD 


<6.0 


<0.4 


<0.3 


<8.3 


2.0 ±0.2 



material selection can significantly reduce the presence of 60 Co in R11410 MOD units, with 
a goal of <0.2 mBq/PMT for complete subdominance. R11410 MOD counting results are 
summarized alongside R8778 results in table [I], and the screening spectrum is overlaid with 
the R8778 spectrum and background in figure |3|. 

The number of PMTs used in an experiment is inversely proportional to the photo- 
cathode area per PMT. Thus a proper comparison of PMT activity must also normalize 
radioactivity by photocathode area. Screening results normalized per unit photocathode 
area are listed for R8778 and R11410 MOD PMTs, as well as for 17 batches of R8520 PMTs 
counted at the University of Zurich and Laboratori Nazionali del Gran Sasso (LNGS) for 
the XENON collaboration |3J], in table @. The R11410 MOD offers by far the lowest activ- 
ity per unit area, with upper limit activities a factor ofXs/x^/XjQ that of the R8520 
positive measurements in 238 U / 232 Th / 40 K, respectively. Activity per unit area for 60 Co 
is measured to be a factor of x ^ below that of the R8520, and would drop dramatically 
after material selection for future R11410 MOD units. 

Table 2. Counting results from table 1, renormalized by PMT photocathode area. Counting results 
from facilities at University of Zurich and LNGS for 17 batches of R8520 PMTs, described in p3 ]. 
are included in summary for comparison. The R11410 MOD offers the lowest radioactivity per 
PMT per unit photocathode area for all benchmark isotopes. In particular, ratios of R11410 MOD 
upper limits to R8520 activities arc xl/3 238 U / xl/4 232 Th / xl/10 40 K. The ratio for 60 Co, 
detected in both PMTs, is xl/3, and can be significantly reduced with further material screening. 
The very low normalized ratios make the R11410 MOD the top candidate for use in the large-area 
tonne-scale LZ detectors. 



PMT 


Phoi 

238 U ( 234m Pa) 


;ocathode-Nor 

238 U ( 226 Ra) 


malized Activit 

232 Th ( 228 Ra) 


y [mBq/ci 

40 K 


n 2 ] 

60 Co 


R11410 MOD 
R8520 
R8778 


<0.19 
<3.4 
<1.4 


<0.013 
0.038±0.002 
0.59±0.038 


<0.009 
0.035±0.002 
0.17±0.019 


<0.26 
2.3±0.07 
4.1±0.13 


0.063±0.0063 
0.18±0.005 
0.16±0.0063 



2.3 Counting Methodology 

Positive signatures of radioactive isotopes are identified through two high-branching ratio 
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energy lines, when available. For measurement of the U chain, lines at 353 keV and 
609 keV are used, corresponding to the decay of 214 Pb and 214 Bi respectively. The 232 Th 
chain is identified by lines at 511 keV and 583 keV from 208 T1 decay. 60 Co is detected 
through peaks at 1173 keV and 1333 keV. 40 K offers only a single gamma line, at 1461 keV. 

It has been noted that 238 U decay chain secular equilibrium can be broken during 
certain manufacturing processes due to the solubility of 226 Ra in water. For this reason, 
a more direct estimate of 238 U concentration is made from measurement of the 1001 keV 
line from 234m Pa, several steps above 226 Ra in the decay chain. The branching ratio for the 
1001 keV line is a factor of x ^ that of the branching ratio for the 353 keV line, and thus 
produces correspondingly weaker upper limits. It should be noted that there is no a priori 
reason to assume a break in equilibrium, unless a positive identification of the 1001 keV 
line is made. It should also be noted that the most relevant, high branching ratio gammas 
are generated in the lower part of the 238 U chain, and this measurement is of the most 
relevance for the estimation of electromagnetic backgrounds. 

Upper limit estimates are always given at 90% confidence level, the standard used 
during LUX material counting. Quoted errors are statistical. Systematic errors are con- 
tributed primarily by uncertainty in sample placement within the chamber and in sample 
self-shielding, and are conservatively estimated to be at the level of ±20%. 

Upper limits are calculated assuming a Poisson-distributed signal. The measured value 
N, in counts, is defined as N = S + B, where S is a positive signal and B is the separately- 
measured background, and either N or B is scaled down such that the number of live days 
for the two quantities is equal. As S depends on the efficiency of gamma detection, it is 
defined as S = £01, where oc is the source radioactivity in units of mBq and £ is the detection 
efficiency already discussed in Sec. |2.1| . £ is measured in counts/mBq, and is obtained from 
Monte Carlo calibrations, detector mass, and counting livetime. For a given confidence 
level CL (equal to 0.9 for all SOLO results), the following integral is numerically evaluated 
to determine the parameter oc': 

CL= J"' p{N,sa+B)da 
Jo p(N,ea + B)da 

where p(N,X) is the Poisson distribution function for measured value ./V and expectation 
A. The value of a' then represents the upper limit source activity at confidence level CL. 

3. Impact on Dark Matter Backgrounds 
3.1 LUX Backgrounds 

Counting results indicating a change of x ^ in 238 U and x g in 232 Th between R8778 PMTs 
and R11410 MOD PMTs are of particular importance, as these isotopes are primary con- 
tributors of neutron backgrounds in dark matter detectors. 2 U and Th both create 
neutrons through (ot,n) reactions in the surrounding PMT materials. 238 U can also un- 
dergo spontaneous fission, directly producing neutrons with a contribution similar to that 
of (o:,n). Neutron rates from (G!,n) are calculated using the Neutron Yield Tool developed 
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Figure 3. SOLO counting spectra for R8778 (14 live days, blue) and R11410 MOD (19 live days, 
red) PMTs, superimposed with a sample background run (21 live days, black). Strong lines in the 
R8778 spectrum indicate the presence of 238 U, 232 Th, 40 K and 60 Co. The reduced activity of the 
R11410 MOD is readily apparent from the lack of distinct features in comparison with the R8778. 



by LUX collaborators at the University of South Dakota [Il5lj . based upon work detailed in 
[llblj , The 238 U concentration within the PMT is assumed to be equally distributed amongst 
all PMT components for the calculations. Neutron yields are calculated separately for each 
chemical compound in each component in order to more accurately account for material 
stopping power effects. 

The (ot,n) rates for the R8778 PMTs, summed for all components, are calculated 
as 0.60 n/PMT/yr and 0.21 n/PMT/yr for 238 U and 232 Th, respectively. The sponta- 
neous fission rate of 238 U is calculated as 1.13 x 10~ 6 n/decay [ li~7f |. yielding an additional 
0.33 n/PMT/yr. This totals to 1.2 n/PMT/yr. The R11410 MOD counting upper limits 
represent a neutron emission rate change of Xyg, or an upper limit of 0.06 n/PMT/yr. 

The LUX experiment could potentially benefit from the use of 62 R11410 MOD PMTs, 
replacing the 122 R8778 PMTs currently in use. Neutron backgrounds are reduced in direct 
proportion to the product of the number of neutrons produced per PMT per unit time and 
the number of PMTs used in the experiment. In the case of LUX, the neutron background 
reduction factor after the R8778 PMTs are replaced by R11410 MOD PMTs is then x^. 
This represents a total of 3.0 neutrons produced from the PMT arrays per 300 live day 
running time; after single-scatter and energy cuts, which reduce event rates by over four 
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orders of magnitude, the neutron background expectation in the WIMP search energy 
region is virtually zero. 

Estimates of electron recoil (ER) background contributions from LUX PMTs have 
been generated using the LUXSim Monte Carlo simulation package developed by LUX 
collaborators [1181] . LUXSim results indicate that the use of the R11410 MOD PMT arrays, 
conservatively assuming activity levels equal to their measured 90% upper limit values, 
would change LUX ER backgrounds by a factor of x-^. This assumes a WIMP search 
nuclear recoil (NR) energy range of 5-25 keVnr 5 matching that used in the LUX proposal. 
It is important to note that this estimate uses the current measurement of R11410 MOD 
60 Co, which is expected to be significantly reduced in future production models. With a 
reduction of 60 Co to levels of <0.2 mBq/PMT, the background reduction factor would be- 
come x^, again making the conservative assumption of 238 U / 232 Th / 40 K levels equal to 
their measured 90% upper limits. The strong reduction in both ER and NR backgrounds 
indicates that the R11410 MOD PMTs can potentially reduce the overall LUX PMT back- 
ground contribution below that of other construction materials used in the experiment, 
creating an overall change in predicted LUX backgrounds of nearly x 4j . 




Figure 4. (Left) LZ ER backgrounds from neutrino scattering and 136 Xe two-neutrino double-beta 
decay, assuming a 2.1 x 10 21 yr lifetime as recently reported in PJ. PMT 238 U and 232 Th ER 
background projections are overlaid for the case of both R8778 and R11410 MOD PMTs. WIMP 
signatures for 100 GeV particles with several possible interaction cross-sections are shown after 
conversion from nuclear recoil to electron recoil energy scales. A 99.5% rejection factor is applied 
for neutrino and gamma spectra. (Right) LZ neutrino coherent scattering backgrounds adapted 
from [liyij . Overlaid are neutron background spectra for both (oc,n) and fission contributions from 
the PMTs. 50% NR acceptance is assumed. 



3.2 LZ Backgrounds 

It is interesting to note that no further reduction in PMT radioactivity beyond the R11410 MOD 
is useful for the LZ experiment. In particular, LZ faces a background "floor" from neutrino 
scattering, causing ER events as well as NR events from coherent neutrino scattering [ ITU ]. 
The spectra for both ER and NR backgrounds are shown in figure |j. Conversion from the 
NR to ER energy scale is performed assuming keV nr /keV ee =0.3 and 50% NR acceptance. 
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ER backgrounds are predicted to be dominated by scattering from solar neutrinos from 
the p-p chain, even with the use of PMTs equivalent to the LUX R8778. With the use of 
R11410 MOD PMTs in the tonne-scale experiment, the PMT NR background contribution 
can be reduced below the level of coherent neutrino scattering from atmospheric and diffuse 
supernova background neutrinos above 5 keV. This will facilitate a direct measurement of 
these scattering spectra. 

4. Conclusion 

The Hamamatsu R8778 and R11410 MOD PMTs are both well suited for use in liquid 
xenon scintillation detectors. The radioactive isotope content of both PMT models has 
been assessed. Counting results for the R11410 MOD PMT yield a combined 238 U + 232 Th 
measurement below 0.7 mBq/PMT. 

Background simulations have been used to estimate the potential reduction in LUX 
backgrounds from replacement of the current LUX PMTs with R11410 MOD PMTs. Sim- 
ulation results indicate that changes of x ^ in LUX PMT electron recoil backgrounds and 
x ^ in neutron recoil backgrounds are expected, conservatively assuming PMT activities 
equal to the 90% upper limits found for all benchmark isotopes, and after reduction of 60 Co 
limits to <0.2 mBq/PMT through material selection in future R11410 MOD units. At this 
point, LUX experiment backgrounds would be dominated by other construction materials, 
and would be well below the level of 0.1 WIMP-like events during the 300 live day running 
time of the experiment. The use of these PMTs in the tonne-scale LZ detector will reduce 
the PMT background contribution below the irreducible contributions from solar and at- 
mospheric neutrinos; no further reductions in radioactivity beyond the measurements for 
the R11410 MOD are necessary for these next-generation experiments. 

The counting results presented here make R11410 MOD PMTs a very competitive 
technology for use in future large-scale xenon detectors. The counting results also offer 
a strong baseline for future photodetector development. The PMTs have a high average 
measured quantum efficiency of 33% at xenon scintillation wavelengths, as well as signal 
gains up to 10 7 at 1500 V. The use of ultra-low background PMTs such as the R11410 MOD 
greatly aids in the development of next-generation detectors with essentially background- 
free discovery potential. 
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